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Abstract 

Explicit solutions for extended objects of a Q-ball type were found analytically in 
a model describing complex scalar field with piecewise parabolic potential in (3+1)- 
and (l+l)-dimensional space-times. Such a potential provides a variety of solutions 
which were thoroughly examined. It was shown that, depending on the values of the 
parameters of the model and according to the known stability criteria, there exist stable 
and unstable solutions. The classical stability of solutions in (l+l)-dimensional space- 
time was examined in the linear approximation and it was shown explicitly that the 
spectrum of linear perturbations around some solutions contains exponentially growing 
modes while it is not so for other solutions. 



1 Introduction 

Q-ball is a nontopological soliton in theories with a global symmetry [TJ. A simple example 
is the model with one complex scalar field (ft with £7(l)-invariant potential V{cft*(ft) which 
satisfies several simple conditions derived in pp. The standard solution for a Q-ball has the 
form 

<ft{t,x) = f{r)e^\ (1) 
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where a; is a constant and f(r) is a monotonically decreasing (in a simplest case) spherically 
symmetric function, tending to zero at spatial infinity. Different applications in cosmology 
(see, for example, [2] for review) encourage investigation of specific features of such classical 
extended objects. However, nonlinearity of classical field equations is a serious obstacle for 
obtaining analytical solutions even in theories with simple polynomial potentials (see, for 
example, |3]), the exceptions like the one presented in [I] (the scalar field potential of this 
model makes it possible to examine analytically even the linear perturbations above the 
Q-ball solution, see [5]) are very rare. Some results of qualitative and numerical analysis of 
Q-ball properties can be found in [6]. 

To obtain an analytic solution the model with parabolic piecewise potential 



V(4>*(p) = M 2 



+ 2ev(v - V(/)*<p)6 ( - 1 

v 



(2) 



was considered in [TJ. Here M 2 > 0, 9 is the Heaviside step function, v is a parameter of the 
model, e is a positive dimensionless constant. Indeed, in both regions f > v and / < v the 
corresponding equation of motion is analytically solvable and the solutions are regular and 
smooth. Matching conditions determine the point r = R such that f(R) = v. Inside the 
sphere with radius R the value of the field / is larger than v and outside this sphere / < v 
and exponentially tends to zero. Thus, the charge and the energy are localized inside the 
ball with the center at the origin r = 0. 

In a solvable model, one can find the whole spectrum of excitations above a Q-ball solu- 
tion, examine the modes responsible for possible instabilities and, in addition, demonstrate 
explicitly the validity of the stability conditions for Q-balls, which can be found in [8| [9| [TO]. 
Indeed, one can examine linear perturbations above a Q-ball solution in analogy with [3], 
where instability of the solution in a theory with potential of the form V = K 2 (f)*(j) — ^(0*0) 2 
was established numerically. The square root in potential (J2]) is an obstacle for analytical 
consideration of linear perturbations above the Q-ball solution. In this paper we propose a 
model describing complex scalar field with potential of the form 

V{<t>*<t>) = M 2 0*0 9 (l - ^f] + (m 2 0> + A)0 - l) , (3) 



2> m 



where M 2 > 0, M 2 > m 2 , 9 is the Heaviside step function with the convention 6(0) 
theories with three and one spatial dimensions. The constant A provides continuity of the 
potential at the point 0*0 = v 2 , A = v 2 (M 2 — m 2 ). This potential is presented in Fig. [TJ Of 
course, this piecewise potential should be considered as a limiting case of some smooth and 
bounded potential. 

Below we will examine Q-ball solutions in a model with potential ((3]) in (3 + 1)- and 
(1 + l)-dimensional space-times and show that the model contains solutions possessing dif- 
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Figure 1: The forms of the scalar field potential described by Eq. ([3]): m? > 0, M/m = 2 
(left plot); m = (middle plot); m 2 < 0, M/\m\ = 2 (right plot). 

ferent properties. As will be shown below, some solutions are stable, some of them are not. 
According to [3], examination of classical stability in the linear approximation will be made 
for the (1 + l)-dimensional case and stability conditions, presented in [HI El EE], will be 
illustrated explicitly. 



2 Q-ball in four-dimensional space-time 



2.1 Analytical Q-ball solution 

First let us consider (3+l)-dimensional space-time. The action of the model has the form 



S 



d x 



du,<i>*d^<f) — M 4>*<j>9 1 



— m 



2 J* 



+ A)9 



- 1 



(4) 



with A = v 2 (M 2 — m 2 ). We consider the standard form of solution ([T]). The simplest regular 
spherically symmetric solution to the equation of motion, coming from fljj), takes the form 



./ 
./ 



2 — 7~i 



B- 



for 



A 



sin ( \/bJ 2 — m 2 7 



for 



f>v< 



(5) 
(6) 



It is clear that if m 2 > 0, then m < \u\ < M; if m 2 = 0, then < \u\ < M; otherwise 
< \u>\ < M. The continuity of the solution and its first derivative at the point r = R such 
that f(R) = v defines the coefficients A and B, which read as 

vR 



B 



.4 



vR 



sin ( y/tu 2 — m 2 R) 



(7) 
(8) 
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and defines the value of the matching radius r = R: 



tan (Vm 2 -m 2 R) = — . (9) 



2 



2 -tu 2 



Since the parameter R is positive, it should be taken to be 



i ( ( 



R = -F= 2 1 arctan - + vr . (10) 

yui 2 — m z \ \ y M z — u z / J 

All the features of our solution can be expressed through this parameter R or, equivalently, 
through the frequency ui. Note that the presented solution is the simplest one - the absolute 
value of the function / is a monotonically decreasing function (which means that the charge 
density is a monotonic function of r as well), which is equal to v only at the one point 
r = R. Obviously, there may exist solutions which cross the line \/0*0 = v several times and 
have nodes (i.e. such that there exist points ^ oo: /(r$) = 0). Because of the simplicity 
of potential (J3J), such solutions can also be found analytically. It would be interesting to 
examine their properties in comparison with those of the simplest solution presented above, 
but this topic lies beyond the scope of the present paper. 

One sees that the scalar field potential is unbounded from below if m 2 < 0. Meanwhile, 
one can always consider a potential which coincides with the one defined by fl3]) for 0*0 < v 2 
but changes its behavior and becomes growing for 0*0 > v 2 , where v > v. In such a case one 
should take only those Q-ball solutions for which the maximum absolute value of the scalar 
field |/(0)| < v. 

Now let us examine the properties of the solution at hand. It is not difficult to calculate 
the total charge and the total energy of the Q-ball. They look as follows: 



« .. , ■ x , , . (M 2 - m 2 ) (Ry/M 2 - u 2 + 1) , 
Q = -i d 3 x( 0*0 - 0*0 = 47T.RW K - - L\^—= =^-L , 11 



) VM 2 



(to 2 — m 2 ) VM 2 — to 2 



E = J d 3 x (^0*0 + ^0*^0 + V(0*0)J = (12) 

, n2 2 2 f (M 2 -m 2 )(RVM 2 -uj 2 + l) \ t An D3 

= A-kR uj v ; H R v (M —m ), 

\ {u 2 - m 2 ) VM 2 - u 2 J 3 v ; ' 

where R is defined by ([TO]) . One can see that E = uiQ + AnR 3 v 2 (M 2 - m 2 )/3 > cuQ. Note 
that the inequality E > uiQ should hold for any Q-ball solution (see simple derivation of this 
fact in Appendix A). Another interesting observation is that 

Q = - V 2 (M 2 - m 2 ) d ^ RZ \ (13) 

CLOU 
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This observation allows one to make another cross-check of our results. Indeed, let us 
differentiate equation (fT2l with respect to the charge Q: 

du „ 47r ,,,o osd(R 3 ) duj 47r 2 2 2 d(R 3 ) did 



dE 
dQ 



= w + ^Lq + — V 2 (M 2 - m 2 )-—t =u + —Q + —v 2 {M 2 - m 2 ' 



dQ" 3 v ' dQ dQ" ' 3 

Using equation ffTBl we easily get the well-known relation [SI El 13] 

dE 



du dQ 



• (14) 



dQ 



UJ. 



(15) 



which should hold for any Q-ball solution. 
2.2 m 2 > 

First, we consider the case m 2 > 0. The E(Q) and E(Q)/Q diagrams for this case for 
M/m = 3 (here and below, without loss of generality, we consider u> > 0) are presented 
in Fig. [2j We see that there are two different branches in this figure, which correspond to 
solutions which are localized in different ways. Localization of solutions with larger energy is 
due to the exponential suppression outside the core [0, R] and we will refer to it as a "wide" 
branch. The size of solutions with smaller energy is just R (see below) and most of their 
charge is concentrated inside the core [0, R]. We will refer to these solutions as a "narrow" 
branch. We see that there is a solution with nonzero minimal charge and energy. This result 
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Figure 2: E(Q) (left plot) and E(Q)/Q (right plot) for the Q-ball in (3+l)-dimensional 
theory, m 2 > 0, M/m = 3. The dashed line corresponds to free particles with E = MQ. 

is very similar to those obtained in [SI [11] where completely different scalar field potentials 
were utilized. Indeed, one can see that the left plot presented in Fig. [2] looks the same as 
the corresponding plots in Fig. 2 and Fig. 3 of [8]. The right plot in Fig. [2] has the same 
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form as the one presented in [JTJ. Of course, these results are also very similar to the results 
obtained in [TJ. 

The quantitative difference between branches can be expressed through the parameter 
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M 2 -uj 2 
(J 2 — m 2 



(16) 



In the case of large Q two different branches, presented in FigfSJ correspond to g ^> 1 and 
g <C 1. For large Q and M ~ m the properties of the solutions are summarized in Table 1. 



Type of solution for large Q 



Asymptotics 



Soliton radius L 



Soliton energy E 



Soliton charge Q 



Wide branch 



9<1 



M 2 —m 2 



v 2 L 



v 2 L/M 



Narrow branch 



0»1 



R 



■Kg 



V M 2 —m 2 



v 2 M A L b 



2 M 3 L 5 



Table 1. Properties of the Q-ball solutions for large Q; m 2 > 0; (3+l)-dimensional 

space-time. 



One can express all values in terms of large Q. For example, if m ~ M the size of wide 
Q-ball is L w ~ and the size of narrow Q-ball is L n ~ ( M f j ) 1 ^ 5 , This difference is the 
origin for the titles of the branches - indeed, for large Q and for M~t)we have L w 3> L n . 
Solutions of the wide branch are very similar to the Q-clouds of [TT] . 

It should be mentioned that the usual thin-wall approximation can not be used for solu- 
tions presented above. Indeed, for the stable branch even in the limit u — > m (in this case 
the solution most rapidly falls off in the region r > R) the solution considerably differs from 
a constant in the region [0, R]. 



2.3 m = 

Now we turn to the second case m = 0. The scalar field potential in this case contains flat 
directions, such type of potentials arise in supersymmetric theories. As can be seen from 
Fig. [31 the properties of the Q-ball solutions in this case look very similar to those discussed 
in the previous subsection. Meanwhile, there is a considerable difference between the cases 
m 2 > and m — 0. In the first case solutions with large charge Q on the lower branch have 
the following energy-charge dependence: E ~ Q. In the second case m = one obtains from 
( JTTT) and ( |T2l) for large Q (i.e for w — > 0) the following energy-charge dependence: E ~ 
The latter relation coincides with the general estimate for potentials of this type, which can 
be found in [2]. 
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Figure 3: E(Q) (left plot) and E(Q)/Q (right plot) for the Q-ball in (3+l)-dimensional 
theory, m — 0. 

2.4 m 2 < 

The third case m? < appears to be completely different from the previous cases. First, 
there are two different phases - solutions in the first phase contain three branches in the 
E(Q) diagram (see Fig. Hj), whereas another phase contains only one branch (see Fig. [5]). 




Figure 4: E(Q) (left plot) and E(Q)/Q (right plot) for the Q-ball in (3+l)-dimensional 
theory, m 2 < 0, M/\m\ = 3. 

The transition between the phases occurs at Mi ~ 0.728. 

V \m\ 

We see that in Fig. H] there exist three branches, two of them intersect each other. There 
also exists nontrivial (i.e. having nonzero energy) solution with the zero charge Q. An 
analogous case was described in [7j for the appropriate range of the parameters of the scalar 
field potential. Note that though the scalar field potential in [7] is bounded from below for 
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Figure 5: E(Q) (left plot) and E(Q)/Q (right plot) for the Q-ball in (3+l)-dimensional 
theory, m 2 < 0, M/\m\ = 0.6. 

any values of the parameters, whereas our potential is unbounded from below for m 2 < 0, 
similar Q-ball solutions having three branches arise in both cases. An interesting observation 
is that the part of E(Q) dependence containing two branches, i.e. the part starting from 
oj = and ending at u = u c < M, where u c corresponds to the solution with minimal 
(this minimum is local in general) energy, resembles the E(Q) dependence found in [5] in 
the model with a completely different potential. The difference between our case and the 
case discussed in [5] is that the E(Q) dependence of [5] has u c — > oo and the corresponding 
branch ends at the point (0; 0) on the Q, E plane. 

The E(Q) dependence of another phase (presented in Fig. [5]) is very similar to the one 
of the model presented in [3]. 

With the help of the explicit solution it can be shown that the thin-wall approximation 
also does not describe the case m 2 < 0. 



2.5 Stability of the solutions 

Now let us discuss the stability of the Q-ball solutions presented above. The first type of 
stability is the quantum mechanical stability. We will focus on the case m 2 > 0. One sees 
from Fig. [2] that the lower branch crosses the line E = MQ, corresponding to free particles 
with the rest mass M, at some charge, say Q x . Thus, for the region of charges Q > Q x the 
energy of the Q-ball is smaller than the energy of free particles and thus such a Q-ball is 
quantum mechanically stable. In the region Q m i n < Q < Q x , where Q = Q m in ~ 2.85^^- for 
M/m = 3, the Q-ball solution is unstable with respect to radiation of free particles. Q-balls 
from the upper "wide" branch are always unstable from this point of view. 
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Now let us consider the Q-ball fission. It is known that Q-balls are stable against fission 
if d 2 E/dQ 2 < (although this condition is known, we present a simple justification of this 
fact in Appendix B). This relation holds for the lowest branches in all three cases (for the 
case m 2 < this is valid for the phase with M > 0.728). Indeed, = |j (see Eq. (ITS])), 
whereas for the lowest branches -t2 < because the charge increases while uj decreases and 
thus d 2 E/dQ 2 < 0. This can also be seen from Fig. [21 Fig. |3] and Fig. HJ 

The last type of stability, which can be discussed here, is the classical stability, i.e. 
stability with respect to small perturbations. The classical stability criterion proposed in 
[HJ [9] implies that a Q-ball solution is classically stable if ^ < (for u > and Q > 0, which 
is exactly our case), a mathematically rigorous proof of this fact can be found in [10] (see 
also references therein). This relation, as it was shown above, holds for the lowest branches 
in all three cases (for the case m 2 < this is valid for the phase with ^ > 0.728) and 
thus we can conclude that Q-balls from these branches are classically stable. An interesting 
observation for the case m 2 < 0, ■My > 0.728 is that classically stable Q-ball solutions lie 
between solutions with locally minimal and locally maximal charges, i.e. there exist stable 
solution with minimal charge (and energy) and stable solution with maximal charge (and 
energy) . 

Meanwhile, it would be interesting to consider linear perturbations explicitly in order 
to examine the upper branches and their possible instability modes. A complete analytical 
analysis of linearized theory in the (3+l)-dimensional case seems to be a rather complicated 
task and it lies beyond the scope of this paper. Nevertheless, we performed analysis of 
perturbations in the linear approximation in a simpler model in (l+l)-dimensional space- 
time for the cases m 2 > and m 2 < 0. The physical properties of the (l+l)-dimensional 
model are very similar to those in the (3+l)-dimensional case and we think that results 
obtained in such a simplified case can be applied to the (3+l)-dimensional well. 



3 Q-ball in two-dimensional space-time 

3.1 Analytical solution and its properties 

Now the action takes the form 

S = j dtdz (d^*d»(j) - M 2 0> e(l- - (m 2 0*0 + A)0 - 1^ (17) 

with A = v 2 (M 2 - m 2 ). A gain we consider the standard ansatz 
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with a dimensionless even function f(z) = /(— z). At the points z = ±R function f(z) is 
equal to v and there are discontinuities in the ordinary differential equation on f(z): 



1 / = 0. 



-d 2 J-u 2 f + M*9 (l-^)f + m 2 9 

Inside the interval (-R, R) the function |/| is larger than v and outside this interval it falls 
off exponentially to zero. The corresponding solution to equation of motion ( TT8|) takes the 
form 



/ 
./ 



ve 



-VM' 2 -u 2 (\z\-R) 



for 



with R defined as 



cos ( y/u 2 — m 2 z ) 
cos(\/u 2 — m 2 R) 

r 1 



for 



f < v 2 , 



\z\ > R 
\z\ < R 



: arctan 



The total charge and the total energy look as follows: 



or 



Q = 2v 2 u 



(M 2 



m 



(Ry/M 2 ^ 



(or 



E 



2v 2 oj 2 



[M 2 - m 2 ) (RVM 




+ 2v 2 (M 2 -m 2 )R, 



(19) 
(20) 

(21) 

(22) 
(23) 



(to 2 - m 2 ) VM 2 - to 2 

where R is defined by ( ]2T]) . Again one sees that E = ooQ + 2t> 2 (M 2 — m 2 )R > ojQ. It 
is straightforward to show that, in analogy with the (3+l)-dimensional case, the following 
relation holds: 

Q = - v 2 (M 2 -m 2 )^^-. (24) 

Using this relation we easily obtain f fT5|) . 

First we consider the simpler case m 2 > 0. The function E(Q) for M/m = 2 is presented 
in Fig.El As in the (3 + l)-dimensional case, there exist two different branches, their 
properties again can be expressed through the parameter g defined by f|T6|) and the effective 
size L. The differences between solutions are summarized in Table 2. 



Type of solution for large Q 



Asymptotics 



Soliton size L 



Soliton energy E 



Soliton charge Q 



Wide branch 



9 < 1 



c\\J M 2 —m 2 



v 2 M 2 L 



v 2 ML 



Narrow branch 



9 > 1 



R 



Kg 



2VW 1 



v 2 M 4 L 3 



! M 3 L 3 



Table 2. Properties of the Q-ball solutions for large Q; m 2 > 0; (l+l)-dimensional 

space-time. 



10 



E 



2mv z 



20 



15 



10 



10 



Q_ 

2v 2 



Figure 6: E(Q) for the Q-ball in (l+l)-dimensional theory (solid line) and for free particles 
with E = MQ (dashed line), m 2 > 0, M/m = 2. 
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Figure 7: E(Q) for the Q-ball in (1+1) -dimensional theory, m 2 < 0, M/\m\ = 2 (left plot) 
and M/\m\ = 1 (right plot). 
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The case m 2 < is also very similar to the one in the (3+l)-dimensional case. Again there 
are two phases, presented in Fig. [71 The transition between the phases in (l+l)-dimensional 
space-time occurs at ~ 1.496. 

3.2 Perturbations and classical stability analysis 

Now let us consider the linearized theory above the Q-ball solution presented in the previous 
subsection. We will be interested mainly in exponentially growing modes which indicate the 
existence of classical instability. We consider small perturbations above the classical solution 
of the form 

4{t,z)=e iut f(z)+e iut Ti(t,z). 
Formally, the linearized equation of motion takes the form 

— rj — 2iufj + u 2 r] + f]" — (M 2 + (m 2 - M 2 )6(f 2 - v 2 )) r]-(m 2 -M 2 )5{f -v 2 )f 2 (rj + rf) = 0, 

(25) 

which can be easily brought to the form 

- t) - 2iuf] + u 2 r] + rj" - (M 2 + (m 2 - M 2 )9(f 2 - v 2 )) r] - FS(\z\ - R) (rj + 77*) = (26) 
where 

. . m 2 - M 2 
F = F(cu) = - 27 

It is worth mentioning that the mixing between 77 and if occurs for potential ()3]) only at 
the points z = ±R through the term with the 5 function and this fact really simplifies the 
consideration. 

There are two obvious solutions to equation (1261) . The first one is the translational mode 

V ~ f{z) (28) 
and the second one corresponds to the existence of the U(l) global symmetry 

7j~if( z ). (29) 

In order to find other possible solutions in the linearized theory we consider the standard 
ansatz for perturbations (see, for example, [31 [5]) 

77 = Vi(^)e i7 ' + ^* 2 (z)e- irt . (30) 

Substituting this decomposition into ( 1261) one obtains the equations 



[-d 2 + U(z)+F5(\z\-R)}^ 1 + F5(\z\- R)^ 2 = (w + 7) 2 ^, 
{-d 2 z +U{z)+F5{\z\-R)}^ 2 + F5{\z\-R)ij l = (u - 7)^2, 



(31) 
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where U(z) = M 2 8(\z\ — R)+m 2 8(R — \z\). It should be noted that formally the case Re7 = 
should be considered separately, because in this case there is no separation of the equation 
( 126]) into terms proportional to e^ Re7 ^ and e~ l ^ Re ^ t , which results in two different equations 
presented above. But it appears that the equation for the spectrum, which will be obtained 
from ( I3ip . is also valid for the case Re7 = 0. 

First, let us consider normalized solutions to equations (l3Tj) . They take the form 

rj( z ) = e 47 V wz a! + e'^e^h, Im/xi > 0, Im/ii > (32) 

with (7 + u) 2 = M 2 + fi 2 and (-7* + u) 2 = M 2 + jj 2 for \z\ > R and 

r,( z ) = e* 7 ' (e^ 2Z a 2 + e"^a 2 ) + e" i7 ** (e r ^ z b 2 + e^b^ (33) 

with (7 + u) 2 = m 2 + fi 2 and (—7* + oj) 2 = m 2 + /t| for \z\ < R. Note that if the conditions 
Im/i! > 0, Im/ix > are not fulfilled, the solutions do not fall off at spatial infinity and thus 
such solutions are not normalizable. 

It is convenient to consider modes which are even and odd in z separately from the very 
beginning. In the first case we have T]'\ z= q = 0. The matching condition at the point z = R, 
coming from ( |26l) . looks as follows: 

v'\z=r+o ~ v'\z=R-o = F(r) + r}*)\ z=R . (34) 

The condition (1341) together with the continuity of rj at the point z = R generate the 
characteristic equation on 7. After some calculations (which are straightforward but quite 
tedious and we do not present the details of calculations here), we can get the following 
characteristic equation for the spectrum: 

%F (G(vi, to)Y(-fe) + G(-jxl -fiDYM) - G(^,^)G(-fLl, -fe) = 0, (35) 

where 

G(/ii, /i 2 ) =112(1- e~^ 2R ) - ^ (1 + e~^ 2R ) , 
Y(ji 2 ) = 1 + e- i2 ^ 2R . 
Note that equation (|35|) can be used only if the conditions 

IH17 ^ 0, (36) 
IH17 = and u - M < Re7 < M - u (37) 

are fulfilled. In this case r] falls off exponentially with z — > ±00 and we have normalized 
solutions. Otherwise (i.e. if Inr/ = and Re7 > M — u> or Re7 < — M + u) one should 
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consider a more general form of perturbations in order to describe modes from the continuous 
spectrum. 

Analogous calculations can be made for the odd modes, for which rj\ x= o = 0. The 
matching condition again has the form (|34p and we get equation (1551) . but now with 

G(ni,fi 2 ) =Hi{l- e~^ R ) - ix 2 (1 + e~ l2 ^ R ) , 

Y(fi 2 ) = -1 + e~ i2 ^ R . 

These formulas also can be used only if (|36|) or (|37|) are fulfilled. 

It is interesting to note that, according to fl36l) and (1371) . the unstable mode (if exists) is 
normalized. As already mentioned, the non-normalized modes from the continuous spectrum 
exist if Im7 = and if M — u < Re7 or Re7 < u — M. 

Now let us turn to an examination of the discrete spectrum of our model. First, it is 
necessary to note that there exist obvious solutions to equation ( 135|) in the case of the odd 
modes. Indeed, if Inry = 0, these solutions to (|35|) are simply ji 2 = or jl 2 = 0, which leads 
to Re7 = — u> — m, Re7 = — u + m, Re7 = u + m, and Re7 = u — m. But it can be shown that 
for these values of 7 the only solution to the initial system of linearized equations of motion 
is 77 = 0. Thus, these roots are unphysical. Note that depending on the value of M/m these 
roots corresponding to unphysical modes may formally lie in the continuous spectrum. 

Equation fl35|) was examined numerically for M/m = 2. We have failed to find solutions 
to (I35I) for Re7 7^ and Inry 7^ 0. For the case IH17 = we have found that solutions to 
equation (1351) may exist (or may not exist) depending on the value of u. The case Re7 = 
and IH17 7^ was examined separately for different values of M/m. We have failed to find any 
exponentially growing odd mode. For u > u c , where u c = u c (M/m) is the frequency of the 
Q-ball solution corresponding to the minimal charge and energy (for M/m = 2 this frequency 
is u c ~ 1.803m), only one exponentially growing even mode was found (it is evident that 
modes with nonzero imaginary part of 7 correspond to classical instability). For uj < u c 
we have failed to find any exponentially growing even mode. Thus, the numerical analysis 
shows that the upper ("wide") branch in Fig. [6] is always classically unstable, whereas the 
lower ("narrow") branch, for which < 0, is always classically stable, which coincides with 
the classical stability criterion of [8j [9]. An example of a nontrivial numerical solution for 
Inry for even excitations together with the function ^2 is presented in Fig. |HJ The results 
of our analysis also show that Q-ball solutions with minimal energy are always classically 
stable. 

The numerical search for exponentially growing modes was also made for the case m 2 < 0. 
We restricted ourselves only to examining the even modes with Re7 = 0. The results of the 
numerical search for Inr/ for even excitations together with the function -Q are presented in 
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-2- 
-4- 

Figure 8: The value of |Im7| as a function of u for the case m 2 > 0, M/m = 2. Solid line - 
numerical calculations; dotted line - the function £\F(u>)\/m. The dashed line corresponds 
to the function 
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Figure 9: The value of |Inry| as a function of u for the case m 2 < 0, M/\m\ = 2. Solid line - 
numerical calculations; dotted line - the function £|F(u;)|/|m|. The dashed line corresponds 
to the function 

aw 

10r 




Figure 10: The value of llnryl as a function of lu for the case m 2 < 0, M/\m\ = 1. Solid line - 
numerical calculations; dotted line - the function ^\F(u)\/\m\. The dashed line corresponds 
to the function 
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Fig. [9] and Fig. [TUJ We see that in the regions where -p- < exponentially growing modes 



are absent, whereas such modes exist if ^2 > 0. The latter means that Q-ball solutions 
in the phase presented on the right plot in Fig. [7] are always classically unstable, because 
= ^ > in this case. This result coincides with the results of numerical analysis in the 
model of [3] in (3 + 1) -dimensional space-time with an analogous E(Q) dependence. 

In the limit of large charges this instability can be obtained analytically in all three cases 
(m 2 > 0, m = 0, m 2 < 0). Let us show it explicitly. Suppose that 7 is purely imaginary and 
I7I ^> M. In this case equation ( |26l) can be rewritten as 

- \l\ 2 V + V" - F5(\z\ -R)( v + 77*) « 0, (38) 

which suggests that the function r\ can be chosen to be real. The even solution to equation 
(138]) takes the form 

rj( z ) = ae l7l *cosh(|7|^), \z\ < R, (39) 
Tj(z) = peWe-M, \z\ > R. (40) 

The continuity of 77 at the points z = ±R and discontinuity of its first derivative at these 
points (see equation (134"]) ) results in 

| 7 | (th{h\R) + 1) = -2F. (41) 

Let us parametrize I7I as 

| 7 | = -£F (42) 
For the solutions at large Q, which are supposed to be unstable, we have oj ~ M and 



„ VM 2 - to 2 
M — m z 

Substituting (021) and (jSj) into (jUJ) we arrive at 

£ (th(£/2) + 1) = 2, (44) 

which results in ^ w 1.2785. The function £|F(a;)|/|m| together with the numerical solution 
for llnryl is presented in Fig. [8] Fig. [9] and Fig. [lOj 

An important remark concerning the validity of the linear approximation is in order. We 
examine classical perturbations, and if the scalar field potential is smooth, we can consider 
linear approximation without any restrictions (we can always choose small enough amplitude 
of perturbations which does not destroy linear approximation). Our case is rather nonstan- 
dard because of the Heaviside step function in the scalar field potential. The 9 function can 
be easily regularized as 



flf^-ll-^ll + th 



v 2 J 2 



v 



(45) 
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where a is a large dimensionless parameter. In the linear approximation it can be expanded 
as 



1 + th 



a 



- 1 



1 



1 + th 



a 



1 



+ 



af(r) + r)*) 



2v 2 cosh 2 (a (j^ 



- 1 



• (46) 



Let us take perturbations at the points \z\ = R. It is clear that the relation BZLtZ? ) <g; i 
should hold in order not to break down the linear approximation. In the limit a — > oo 
(which is exactly the case for which the linearized equation (1261) was obtained), this relation 
obviously leads to the constraint rj + 77*|ui=r = 0, otherwise formally we get the breakdown 
of the linear approximation. Thus, let us consider the linearized theory with the constraint 



r] + r]* | \ g \ =R = 0. 



(47) 



Substituting (1301) into (JITl) and taking into account the fact that the relation ff4T|) should 
hold at any moment of time, we get 

ifa{±R) = -M±R)- (48) 
Now we take the first equation of (13T1) . Taking into account ( l4~8l) . it can be rewritten as 

U> x = (co + 7)^ (49) 



and 

where the operator L is defined as L = —d 2 z +U(z). Multiplying 
the coordinate z and using ( 1301) . we get 

(c + 7 ) 2 = (u; + 7 *) 2 . 



(50) 

9j by if)* , integrating over 



(51) 



The latter equation has two solutions. The first one leads to Im 7 = 0, which means that 
exponentially growing modes, indicating instability, are absent in this case. The second 
solution is 7 = — to + iji, where 7i = Inry. Using (T49l) . we easily obtain 



(Vil^lVi) = -T^il^i)- 



(52) 



But this relation can not be fulfilled, because (if>i\L\if>x) is nonnegative for any if)\. Indeed, 
the eigenfunction of the operator L, corresponding to the eigenvalue u 2 > 0, is f(z). This 
eigenfunction has no nodes, which means that it is the lowest mode and thus all other 
eigenvalues are also larger than zero, leading to non-negativity of (ipi\L\ipi) . 



18 



Analogous calculations, leading to the same result, can be performed for the function ^>%. 
Thus, we have shown that formally there are no modes indicating instability if constraint 
f )47|) is imposed. Meanwhile, this constraint is a consequence of the existence of a generalized 
function in our potential. A physically reasonable potential should be smooth, leading to 
nonzero, though sometimes very small, amplitude of perturbations. It this case in the lin- 
earized equations of motion (126]) delta-function should be replaced by some smooth function 
containing some parameter of regularization a (like the one in fj45l)). For very large, but 
finite a, the linearized equation of motion, as well as corresponding solutions for perturba- 
tions and eigenvalues 7, may look very similar to those in the a — > 00 limiting case. Of 
course, most probably such an equation can not be solved analytically. But the difference 
between the exact solution and the solution in the a — > 00 limiting case is controlled by the 
parameter ~ and for very large a the corrections to the case a — > 00 are supposed to be 
very small. This reasoning justifies the use of the "relaxed" linearized theory described by 
equation ( 126|) . without constraint f|47|) . This "relaxed" theory allows one to see what could 
happen with the stability in a more realistic case of a smooth scalar field potential. 

4 Conclusion 

In the present paper we discussed Q-ball solutions in a model describing complex scalar field 
with a piecewise parabolic potential ((3]). It was shown that due to the simplicity of the 
potential Q-ball solutions can be found analytically, which really simplifies the analysis of 
their properties - the spectra of solutions were obtained analytically in (3+1)- and (1 + 1)- 
dimensional space-times for different values of the parameters of the model. For uj — > M the 
Q-ball solutions are very close to the condensate line E = MQ and this fact can be interesting 
for examining the Q-ball production. It should be noted that the Q-ball solutions presented 
in this paper can not be described by the standard thin-wall approximation. 

The stability of the obtained solutions was also examined. In the simpler theory in 
(1 + l)-dimensional space-time it was shown explicitly that solutions with ^2 > are always 
unstable, whereas solutions ^ < do not contain exponentially growing modes, at least of 
the simplest form ~ e' 7 '*. 

Though the potential ([3]) is very simple, it provides a variety of Q-ball solutions of 
different types. We hope that the existence of exact and simple analytical Q-ball solutions 
with known properties (such as their stability and E(Q) dependence) allows one to consider 
the model, presented in this paper, as a useful tool for examining different phenomenological 
scenarios involving Q-balls. 
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5 Appendix A 

Let us prove that inequality E > ujQ holds for any Q-ball solution. Consider (D+l)- 
dimensional space-time. According to (JTJ we can use the time-independent effective action 

Seff = J d D x (u 2 f 2 - dtfdif - V(f)) , (53) 

where i = 1, D, instead of the initial one. Suppose that there exists a solution /(x) to the 
corresponding equation of motion. Let us apply the scale transformation /(x) — > f\(x) = 
/(Ax) to this solution and substitute f\(x) into the effective action instead of /(x) (this 
technique was used in [12] to show the absence of time-independent soliton solutions in some 
models with a nonlinear scalar field). Now it takes the form 

*e// = / d D x ^f(Xx) - A 2 ^y/(A^)^y/(Ax) - W(Ax))) (54) 



J ^d D x {u 2 f - X 2 dJdJ - V(f)) 



where we have passed to the new variables Ax — > x in the second integral. According to the 
principle of least action we have 

S x 

eff - (55) 



dX 

and thus 



A=l 



J d D xV(f) = u 2 J d D xf 2 - J d D xd l fd l f. (56) 

Substituting the latter equation into the definition of energy and taking into account the 
definition of charge Q = 2u J d D xf 2 , we get 

E = J d D x {u 2 f 2 + difdif + V(f)) =ujQ + ^J d D xdJdJ > coQ. (57) 
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6 Appendix B 



Let us show that if the condition d 2 E/dQ 2 < is fulfilled, then the Q-ball is stable against 
fission, i.e. against a decay into Q-balls with smaller charges. Suppose E(Q), where Q is 
supposed to be nonnegative, is a positive monotonically increasing function in the region of 
charges we are interested in. In this case condition d 2 E/dQ 2 < means that 



E(Q) 



dQ 



E(Q) 



dQ 



(58) 

Q=Qi+Q 



for Q > 0. Let us integrate this inequality over the coordinate Qi in the region [Qi, Q 2 ]. We 
get 

E(Q 2 ) - E(Q 1 ) > E(Q 2 + Q) - E(Q 1 + Q), (59) 
which can be rewritten as 

E{Q 2 ) + E{Q 1 + Q)> E{Q 2 + Q) + E(Q X ), (60) 

If we take Q 1 = and if E(0)=0, then (JSDJ) leads to 

E(Q 2 ) + E(Q)>E(Q + Q 2 ), (61) 

which means that Q-ball fission is energetically forbidden. But in many models including 
the one discussed in the present paper there exists a minimal charge Q m in ^ such that 
E(Qmi n ) = E min 7^ 0. In such a case, one can try to redefine the function E(Q) in the region 
[0, Qmin) i n order to get an auxiliary function E aux {Q): E aux (0) = 0, E aux {Q) is a smooth 
monotonically increasing function for Q > 0, d 2 E aux {Q)/dQ 2 < and E aux (Q) = E(Q) for 
Q > Qmin- If it is possible to construct the function E aux (Q), then inequality ( 1oT|) is valid 
for Q, Q 2 > and, consequently, for Q, Q 2 > Q m i n (of course, any Q-ball with the charge 
Q < IQmin is always stable against fission regardless of the sign of d 2 E/dQ 2 ). 

To show that one can always construct such an auxiliary function for a Q-ball let us 
consider the function E(Q) = uQ, where the constant oj is defined by u = E min /Q r 



'mm- 



Recall that ^M\ n _ n = co m in- It is evident that if §IL_ n = u > u min = ^\ n _ n then 
one can always construct a function E aux (Q), otherwise it is impossible, see examples in 
Fig. EH 

Using equation ([57]) we get 

2 r 

and thus Co > u m i n - The latter means that we can always construct a function E aux (Q) and 
inequality (JSQ fulfills if d 2 E/dQ 2 < is valid. 
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Figure 11: Solid lines correspond to functions E(Q) and E(Q), dashed line corresponds to 
continuation of the function E(Q) (i.e., to the function E aux (Q)). 
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